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Abstract 

Study of parton distribution function (PDF) is a topic of significant 
interest in QCD. To overcome the shortcomings of conventional PDFs, 
several alternative methods have been suggested in recent years. The 
present work reports the x-distribution of the non-singlet structure 
function using the complete solution of Taylor-approximated DGLAP 
equation. The formalism is more general than the similar ones re- 
ported in recent literature. Our predictions are compared with the 
CCFR neutrino data as well as exact results. 
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1 Introduction: 



Recently pQ we have reported the t-evolution of the non-singlet structure 
functions at low and high x as the complete solution of Taylor-approximated 
DGLAP equation [2J [21 HI [5, 6J. The present paper reports the prediction of 
corresponding x-evolution. Although the conventional wisdom on DGLAP 
evolution does not favour such a feature, such a possibility was reported 
[7] sometime back. The main reason is that in this approach, the DGLAP 
equation is converted into a first order differential equation in two variables 
x and t ( t = log instead of only t and has x evolution as its natural 
outcome. 

There are, of course, various methods that have been developed recently 
like the brute-force [HIE], Laguerre-polynomial [TO], [XT] and Mellin-transformation 
[T2l [T3l [HI [THl [16] etc. for numerical solution of DGLAP equation. All these 
methods have good numerical accuracy and computational efficiency, par- 
ticularly the Mellin-transformation one. All these solutions are very useful 
in the studies of a top-down scenario in order to determine the origin of 
ultra-high energy cosmic rays. In addition, the matrix technique to get the 
numerical solution of the DGLAP evolution equation has also been proposed 
|17j . which provides very good precision. 

In recent years, the study of parton distribution function (PDF) has it- 
self become a topic of significant interest. The standard approach to the 
problem is based on a choice of a specific functional form in x with a few 
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parameters and to obtain fit of them by comparing the exact solution of the 
evolution equation with global data. Based on this approach, several sets of 
pdf have been produced in recent years [TBI [HI EQ1 [2~H More recently, 
however, shortcoming of such an approach has been pointed out : the form 
of parametrization is itself a source of bias. To overcome this limitation, nu- 
merical method incorporating the notion of neural network [23"l 12^1 121)112(31 I2T] 
has been suggested and pursued. Alternative conventional models like sta- 
tistical [281 [221 [30] and chiral quark model [3TJ have also received cosiderable 
attention in this respect. 

The present paper reports the analytical ^-evolution of the non-singlet 
structure function in the formulation of Ref. [Tj which is much more general 
than that of earlier works [7] and compares with the other approaches. Since 
such an analytical solution is in good agreement with experimental data, 
hence this type of study bears important significance. In Sec £j2j we develop 
the essential formalism while in Sec $31 we will discuss the results. Sec §H 
will summarize a conclusion. 



2.1 Particular Solution of DGLAP equation at low x: 

The DGLAP evolution equation for non-singlet structure function in stan- 
dard form is given by [5] 



2 FORMALISM: 



8F NS (x, Q 2 ) 
d log Q 2 



Pns (x,Q 2 )®F ns (x,Q 2 ) 



(1) 
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where x is the usual Bjorken variable 

and Q 2 is the four momentum transfer in deep inelastic scattering [32] pro- 
cess. The symbol (g> stands for the usual Mellin convolution in the first 
variable defined as 



/>'(;)■ 



a(x) (g> f(x) 

The non-singlet kernel Pns{x,Q 2 ) has perturbative expansion in orders of 
the strong coupling constant a s (Q 2 ) as 



P Ns{x , q>) = ^i#> (l ) + f ^) 2 Pg (l ) + (4) 



where P^ s (x) and P^ s (x) are the one loop [5] and two loop [33], EH [35], [36] 
contributions to the non-singlet structure function respectively. In case of 
LO, we neglect two loop contribution and so eq. ([1]) can be rewritten as 



Since the evolution of the non-singlet structure function is independent of 
the gluon, hence one can write from earlier equation, 

^.f)p, ?(l)8 ^ (l , c1 (6) 

where we have [15 
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1 - x + -5(1 - x) 



(l-x)+ 2 



(7) 



where Cp is given by Cp = 2j^~ an d * s defined by 



dx 



/(*) 



(l-x) 



/(*)-/(!) 



(8) 



Here f(x) is an arbitrary function and Nc is the number of colours. 

Using this explicit form of the splitting function Pj '(i), eq. flSj) can be 
written as 1381 



dF NS (x,t) a,(t) 



dt 



2tt 



-{3 + 41og(l-x)}F^(M) 



(9) 



Here t = log (jj) an d the running coupling constant in Leading order (LO) 



is a s (t) -- 
with 



An 

Pot 



33— 27V 

— 2 — iV/ being the number of quark flavours. 
This expression can be re-written as 



8F Nb (x,t) A f 



dt 



{3 + 41og(l- x)}F NS (x,t) 



+2/^ {(1 + z 2 )F NS (Z,t\ - 2F NS (x,t)\ 



(10) 



nere /ly - 3ft) - 33 _ 2Nf ana ^ — t ■ 



For small x (x -C 1), it is justified to neglect higher derivatives of aF 
so that F NS (j, tj occured in the R.H.S. of eq. ffTQ]) can be written as 

k dF NS (x,t) 



F NS (~,t)=F NS (x,t) + xJ2u" 



k=l 



dx 



(11) 



where u = 1 — z. 

Putting eq. ffTTj) in eq. (jTUI) and performing the u-integration yields the 
Taylor-approximated DGLAP equation as 

Q(x, t) 9 ^^ + P(x, t) dFNS d ( x X,t) = R (*> *) ( 12 ) 

where P(x,t), Q(x,t) and R(x,t) are explicitly mentioned in of Ref. [1] as 



-A f x 



21og(±) + (l-x 2 )] Pl ( x ) 



P{x, t) = —± ^ i = (13) 



g(x,t) = i (i4) 

and 

R{x t) = A / [3 + 41og(l-x) + (x-l)(x + 3)] = (15) 

Eq. ffT2l is frequently referred to as Lagrange's equation [32J HO]. Its 
general solution is obtained by solving the following auxiliary systems of 
ordinary differential equations 

dx dt dF NS (x,t) 



P(x,t) Q(x,t) R(x,t)F NS (x,t) 
which subsequently gives 

dx dt dF NS (x,t) 



(16) 



(17) 



A Or) Q l (t) R l {x)F NS {x,t) 

where P\{x) = tP(x,t), Qi(t) = tQ(t) and Ri(x) = tR(x,t). The most 
general form of the solution of the Taylor-approximated DGLAP equation 
(eq. (H2D) is [I], 

v = au n{x ^ + (3 (18) 
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where n(x, t) is any real function of x and t and j3, a are two arbitrary 
constants, u and v are two explicit functions of x and t given by [I], 



= tX NS (x) 



and 



(19) 



(20) 



= F NS (x, t)Y NS (x) 
which are two independant solutions of eq. (TTT1) . 

In order to develop our formalism, we take the numerical value of n(x, t) = 
n num = n and use the following two physically plausible boundary conditions 



F NS (x,t) = F NS (x ,t) (21) 

for some low t = to and 

F NS (l,t) = (22) 

Putting the boundary conditions (12T1) and ( 1221) in eq. (JT51) . along with ex- 
pressions for u(x,t) and i>(x,t) from eq. (fT9l) and eq. (!20|) respectively, one 
gets 

F^x, *)Y* 5 (x) = at n {[X NS (x)] n - [X NS (l)] n } (23) 

Defining 



F NS (x ,t)Y™(x ) = at n {[X N »(x )] n - [X» s (!)]»} (24) 



-TVS 



where X (x) and Y (x) are given as [T] 



X^x) = exp 



dx 



(25) 
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and 

Y NS (x) = exp 

This results in 

F NS (x,t) = F NS (x ,t) 
As we have pQ 



R^x) 



dx 



(26) 



Y NS (x ) [X NS {x)] n - [X NS {l)] n 
Y NS (x) [X NS (x )} n - [X NS (l)] n 



X Nb (l) « 0, 



(27) 



(28) 



hence eq. ( 12T1) yields for n > 0, 

F w5 (x,t) = J P w5 (x ,t) 



r^Cxo) [x NS (x)] n 



(29) 



F ws (x) [X W5 (x )] n 
On the other hand for n < 0, the R.H.S of eq. (J27J) will involve inverse 
of X^^l) which is singular and ill-defined, and hence excluded on physical 
grounds. Putting the values of X NS (x) and Y (x), we get x distribution 
obtained from complete solution at low x, 

-n Ri(x) 



F" b (x,t) = F Nb (x ,t) exp 



P 1 (z) + K(x)| dx 



(n > 0) (30) 



where P\{x) and R\{x) are defined as [T] 



P^x) = -A/a; 



2 log 



x 



'l-x 1 



R x {x) = A/[3 + 41og(l - x) + (x - l)(x + 3)] 
Eq. ( 1301) is also true for quark of each flavour individually i.e 

— n R±(x) 



(31) 
(32) 



lim qAx, t) = qi(x , t) exp 

x^0 



■CO 



P 1 (x) + P 1 (x) \ dX 



n > 0) (33) 
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Combining t-evolution for quark PQ at low x, viz, 



q(x,t) = q(x,t ) 



to 



we obtain 



linages,*) = qi(x ,t ) I — J exp 

x^O \t ( 



(n>0), 



-n -Ri(x) 



(34) 



Pi(x) P x (x) 



(n > 0) (35) 



containing t-evolution as well. Eq. (l30l) and eq. (1351) are more general than 
the similar x distribution reported recently 



2.2 Particular Solution of DGLAP equation at high x: 

For high x, the most general form of the solution of the Taylor-approximated 
DGLAP equation is pQ 



v = au +fj 



(36) 



where p(x,t) is another function of x and t similar to n(x,t) in eq. (fT8l) . As 
in case of low x, we assume p(x, t) = p num = p so as to yield 



F NS (x,t) = F^(x ,t) — 



NS, 



Y'"°{x ) [X'"°{x)] p - [X' ns {1)]p 



Y' Nb (x) [X' Nb (x )]P- [X' ns (1)]p 



(37) 



where X' NS (x) and Y'" a (x) are defined as pQ 



NS 



and 



X' NS (x) 



Y' NS (x) 



exp 



dx 
P[(x) 



exp 



Ri(x) 
P[(x) 



dx 



(38) 



(39) 
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As X >NS (1) is singular [Ij, for p > 0, R.H.S of eq. (1571) is ill-defined and 
excluded on physical grounds. If on the other hand, for p < 0, R.H.S of eq. 
( 157)) involves X /JVS (1) , which vanishes identically. Under this condition 

(p<0), 

N S N S 

Putting the values of X' (x) and Y' (x), we finally get the x-distribution 
obtained from complete solution at high x to be 

RAx) 



F NS {x,t) = F NS {x ,t)exp 
where 



-p 



P{(x) P[{x) 



dx 



P[{x) = -A f x(x - l)(x 2 + x + 4) 
3 

Eq. ()4ip is also true for each quark flavour individually. 



lim qi(x,t) = qi(x Q ,t) exp 

(1— a:)— +0 



~ p +^M\dx 



X'() 



P[{x) 



Combining t-evolution for quark pQ at high x, viz, 



(p < 0) (41) 



(42) 



(P < 0) (43) 



t) = g(s, t ) (^— J (p<0), 

we obtain 



(44) 



lim ftfot) = qi(x ,t ) — exp 

(l-x)->0 \Iq, 



t \ p 



P[{x) + P[(x)) dX 



(p<0) 
(45) 



containing t-evolution as well. 
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2.3 General Solution of DGLAP equation at low and 
high x: 

In Subsection 12.11 and Subsection 12.21 , we considered the numerical values of 
n(x,t) and p(x,t) that differentiated the low x from the high x regimes in 
case of x-evolution. Now we attempt to find out the general and complete 
analytical solution of DGLAP equation using the explicit functional forms 
of n(x, t) and p(x, t) and their x — > 1 and x — > xq limits at low and high x 
respectively. 

Now in case of low x, considering the functional form of n(x, t) and pro- 
ceeding as in Subsection 12. 1[ we get 

F NS (x, t)Y NS {x) = a{[tX NS (x)] n ^ - [tX NS {l)] n ^} (46) 

Defining 

F NS (x ,t)Y NS (x ) = a{[tX NS {x )] nM - [tX NS {l)] n ^} (47) 

we get, 



F Ns (rt) = F Ns (r f) Y NS (*o) [tx^( X )r^-[tx^(i)r^ 

{ , ) 1 °' ' Y NS (x) [tX NS (x )] n (*o,t)-[tX NS (l)]"W K } 

Now proceeding as in Subsection 12.11 we ultimately get the general x- 
distribution at low x, 



F NS (x,t) = F«»(x ,t) 



NSi 



Y NS (x ) t<*$ [X NS {x)] n ^ 
Y NS (x) H*"-*) [X NS (xo)] n ^ ' 
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(n(l,t) >0) (49) 



instead of eq. (1251) . 

Similarly, considering the functional form of p(x, t) and proceeding as in 
Subsection 12.21 we ultimately get, 



Y' NS (x) #»(*<>,*) [X' NS {x q )]pM' 

instead of eq. (HOI . 

Eq. (Pf9l) and eq. (150]) can be re-expressed as 



(p(l,t)<0) (50) 



F NS {x,t) = F JV5 '(x ,t)exp 



-j-n{x,t) 
fn(x ,t) 



X- 



exp 


— n(x, t) j 


f — -—dx 

Pi(x) 


exp 


-n(x,t) j 


f — -—dx 
Pi(x) 



X=XQ 



and 



{n(l,t) >0) 



(51) 



F NS (x,t) = F JVi '(x ,t)exp 



ATS/ 



i?i(a;) 

y 7^) 



exp 


-P(x, t) J 


1 ' Pi (x) ix _ 




exp 


-p(x,t) J 




X=XQ 



(p(l,t) <0) 



(52) 
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We now put forward the generalized statement as: The most general 
solution of the Taylor-approximated DGLAP equation at low and 
high x are obtainable in terms of two arbitrary real functions n(x, t) 
and p(x, t) such that n(l, t) is greater than zero and p(l, t) is less than 
zero for any t. In case, n(x,t) and p(x,t) have only mild variation in 
t, then n(l,t) ~ n(l) > and p(l,t) ~ p(l) < can be considered as 
the x — > and (1 — x) — > limit of an interpolating function H(x,t) 
to be determined empirically from data. 

2.4 Alternate Solution of DGLAP equation at low and 
high x: 

Alternatively, we can also determine x-evolution of F NS (x,t) if a and (3 of 
eq. (fT8|) satisfies additional relation (3 = a m [HJ 01], where m is a positive 
integer. In this case, for low x, we get from eq. ([15]) . 

^(x, t)y^ 5 (x) = a[tX Af5 (x)] n + a m (53) 

where eq. ( |T9l) and eq. ( 120|) have been used. 

Differentiating eq. fl53l) with respect to a, we get 

a = -( — j t—{X NS (x)]— (54) 

which forbids m=l. 

Putting this value of a in eq. (1331 . we ultimately get 

F^M^Or) = j- (1)^ - (J^^Utt[X N3 {x)]& (55) 
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we get, 




(x ,t) 



Y NS (x ) [X NS (x)}^ 



(57) 



If m is very large, -p^r ~ 1, then eq. (|571) reduces to eq. (1291) and equiva- 
lently eq. (J30l) . However, unlike in the previous case, there is no additional 
condition of n > 0. Similarly in case of high x too, we get the ^-distribution 
from the alternative way, as in low x, similar to eq. fT4T)l) and eq. (1411) . but 
without the additional condition on p i.e p < 0. Thus unlike in the general 
case, the alternate solution cannot distinguish between low x and high x 
regimes, which is its limitation. 

3 Results: 

3.1 Result on the basis of Particular solution: 

We obtained an interpolating function H(x, t) from CCFR Neutrino data [1] 
viz., 



such that for x — > and x — > 1, it has desired behaviour of n(x, t) and p(x, t) 
respectively, yielding n num < 0.3586 and p num > —2.151 Assuming these 
limits as process independent, we use them to study x-distribution in Figure 
[]Ja-e) and Figure [2]^a-e) for low and high x respectively using the following 
equations that have been redefined from eq. (129]) and eq. (T40l) respectively 



H(x, t) = 0.382[-5.976x + 0.996(1 - x)} 



(58) 
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as 




Y NS (x ) [X ns (x)} h ^ 
Y NS {x) [X NS (x Q )] H M 



(59) 




Y' NS (x ) [X' ns (x)] h ^ 
Y' NS {x) [X' N3 (x )] H ^,t) 



(60) 



We compare our results with CCFR data [IS, H6] as well as GRV(98) 
03 SB] and MRST LO 2001 [2Ql EH [22] exact results. For getting the exact 
results, we take the non-singlet structure function in leading order as |49j . 

xF 3 (x,t) = -x(u + d) + x(d + u)[\v ud \ 2 + \v cd \ 2 } + 2xs[\v us \ 2 + \v cs \ 2 } - 2xc. 



where u = u(x, t) etc. and v ud , v us etc. are the relevant CKM matrix elements 
[5D]. We have taken v ud « 0.97, v us « 0.22, v cd w 0.22 and v cs w 0.97. 

For low x (Figure [T](a-e)), our preditions are found to be tallying very 
well (Within 3% limit) with both data and exact results in the lower part 
of the x spectrum, but it does not match so well in the upper part of the x 
spectrum. 

For high x (Figure [21^ a-e)), on the other hand, our predictions are closer to 
both data and exact results (Within 5% limit) in the extreme lower part of the 
x spectrum. In the middle part of the x spectrum, coincidence between our 
predictions and both data as well as exact results is not so well as expected. 
For extreme higher parts of the x spectrum, there is excellent coincidence 
between our predictions and both data and exact results (Specially in the 
region 0.60 < x < 0.75). 



(61) 
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The deviations of our predictions from data and exact results is pre- 
sumably due to the neglect of finite x corrections coming from the higher 
derivatives of dF d ^ x,t ^ in the Taylor approximation (Neglected in eq. (3) 
and eq.(35) of Ref. [I]). It may also be partly due to the neglect of explicit x 
and t dependence of n(x, t) and p(x, t) in developing the seperate formalism 
of low and high x. 

3.2 Result on the basis of General solution: 

We redefine eq. f )49l) and eq. fl50|) as 



\x,t) r \x a ,i) Y , nS(x) ( B M [X'" s (i f(-..')' lM ' 
which are valid for low x and high x regions respectively. 
On the basis of the above equations, we study the ^-distribution in Figure 
Ufa-e) and Figure [2](a-e) for low x and high x respectively. The dash line in 
the above set of figures represents the result obtained on the basis of general 
solution. 

For low x, our predictions based on the particular and the general solution 
does not differ appreciably in the lower part of x-spectrum, but the general 
solution is found to match better with both data and exact results in the 
upper part of x spectrum. On the other hand, for high x, the general solution 
matches better than particular solution with both data and exact results in 

16 
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Figure 1: (a-e) xF 3 versus x in bins of Q 2 obtained from eq. fl59|) (Solid line) 
and eq. fl62|) (Dashed curve) for low x values. This is compared with the 
GRV98 H3 SB] (Dotted curve) and MRST2001 [HI E2] (Dashed dotted 
curve) exact results calculated using eq. ([Ml) . This is also compared with 
CCFR data HHIS]. 
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Figure 2: (a-e) xF 3 versus x in bins of Q 2 obtained from eq. (160]) (Solid line) 
and eq. fl63|) (Dashed curve) for high x values. This is compared with the 
GRV98 HUSH] (Dotted curve) and MRST2001 [201 EH E2] (Dashed dotted 
curve) exact results calculated using eq. flBTT) . This is also compared with 
CCFR data 03 HE]. 
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the lower part of x spectrum, but it is found that the opposite is true for 
higher part of x-spectrum. 

Thus we find that acceptance of explicit x and t-dependence of n(x, t) 
and p(x, t) gives better result in the higher x spectrum of low x and lower x 
spectrum of high x. 

3.3 Incorporation of large x effect: 

For low x , we find in Subsection 13 .21 that predictions based on the most gen- 
eral solution have shoot up the CCFR neutrino data as well as exact results 
in the higher part of x -spectrum. This overshooting can be substantially 
reduced by the introduction of large x effect. Such finite x corrections, which 
may be considered as the conglomeration of the terms neglected while getting 
eq. (fTTj) . can effectively be simulated by an empirical multiplicative damping 
factor (1 — x) r , (r is positive) in the structure function (General Solution) 
in eq. (1621) . We have varied r from 4 to 9 and found r around 7 to yield 
the best fit to the data and exact result. This is clearly shown in Figure 
[31 where we have plotted this modified structure function at different x for 
fixed Q 2 = 7.9 GeV 2 . Similarly, we multiply the structure function (partic- 
ular Solution) in eq. fl59l) by multiplicative damping factor (1 — x) r and plot 
this modified structure function versus x for fixed Q 2 = 7.9 GeV 2 in Figure 
HJ This time we find the value of r around 8 gives the best fit to the data, 
clearly pinpointing the supremacy of the general solution over the particular 
solution. It is pertinent to note that the modified structure function damps 
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Figure 3: Modified xF 3 versus x for fixed Q 2 = 7.9 GeV 2 simulated by an 
empirical multiplicative damping factor (1 — x) r in the structure function 
(General Solution) in eq. f|B2|) for low x values. This is compared with the 
MRST2001 PDl |2Il |22] (Dashed dotted curve) exact results calculated using 
eq. (EU) as well as CCFR data 03 H6] . 
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Figure 4: Modified xF 3 versus x for fixed Q 2 = 7.9 GeV 2 simulated by an 
empirical multiplicative damping factor (1 — x) r in the structure function 
(particular Solution) in eq. (|59p for low x values. This is compared with the 
MRST2001 [201 EH E2] (Dashed dotted curve) exact results calculated using 
eq. (EU) as well as CCFR data 03 US] . 
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considerably in the higher part of x -spectrum (of low x) only due to the 
damping factor. In the lower part of the x spectrum (of low x), there is very 
little damping because (1 — x) r — > 1 when x —>■ 0. 

3.4 Comparison with other works: 

Let us now discuss how the present work differs from that of Ref. [7] , [43] and 
[Sj. In Ref. [7j, Taylor-approximated leading order coupled DGLAP equation 
for singlet (F 2 5 (x,t)) and gluon distribution (G(x,t)) were solved using the 
relationship G(x,t) = k(x)F2 (x,t)) and the particular solution v = au + (3 
(a and (3 are constants) implying n(x,t) = 1 in the general solution of the 
present formalism (eq. (Tl8|) ). The results are then compared with deuteron 
data using the relation F^(x,t) = |F 2 5 (x,t). 

In Ref. [43], a similar analysis was done for deuteron structure function 
using the additional relation (3 = a 2 . In Ref. jH], the analysis was further 
extended to Next-to-Leading order for deuteron, proton and neutron struc- 
ture function assuming the relationship (3 = a m (m = 2,3,4,5 ) implying 

that x evolution of the formalism is not unique. In the present work, we 
have made the formalism much more general abandoning the linearity of u 
but still preserving that of v (which is equivalent to structure function) in 
the general solution. 
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4 Comments and conclusions: 



To conclude, we have shown that the x-distribution of the non-singlet struc- 
ture function obtained from Leading Order Taylor-approximated DGLAP 
equations are compatible with data [HJ EJ6J as well as exact results for low 
and high x [201 [221 H7[ HH] in selected parts of the x spectrum. We have 
also shown how the introduction of damping factor considerably brings the 
modified structure function closer to exact resuts in case of low x. We have 
also discussed how the present formalism is an improvement and generaliza- 
tion over those of Ref . [3 M, M, EH] ■ 
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